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ABSTRACT 

Aims. The consequences of a newly suggested value for the SN explosion energy E Bn = 1.2 X 10 51 erg are explored for the case of Tycho's 
supernova remnant (SNR). 

Methods. A nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) is employed to investigate the properties 
of Tycho's SNR and their correspondence to the existing experimental data. 

Results. It is demonstrated that the large mean ratio between the radii of the contact discontinuity and the forward shock is consistent with the 
very effective acceleration of nuclear energetic particles at the forward shock. It is also argued that consistency of the value _E sn = 1.2 x 10 51 erg 
with the gas dynamics, acceleration theory, and the existing 7-ray measurements requires the source distance d to be greater than 3.3 kpc. The 
corresponding ambient gas number density is lower than 0.4 cm" 3 . Since the expected 7-ray flux strongly depends on the source distance, 
F 7 oc d , a future experimental determination of the actual 7-ray flux from Tycho's SNR will make it possible to determine the values of the 
source distance d and of the mean ambient gas density. A simple inverse Compton model without a dominant population of nuclear CRs is not 
compatible with the present upper limit for the 7-ray emission for any reasonable ambient interstellar B-field. 

Conclusions. Given the consistency between acceleration theory and overall, as well as internal, gas dynamics, a future 7-ray detection would 
make the case for nuclear particle acceleration in Tycho's SNR incontrovertible in our view. 

Key words. (ISM:)cosmic rays - acceleration of particles - shock waves - supernovae individual(Tycho's SNR) - radiation mechanisms non- 
thermal - gamma-rays:theory 
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1. Introduction 

Nonlinear kinetic theory of diffusive CR a c celer ation in SNRs 
dBerezhko et all 1 19961 iBerezhko & Volkl 1 1997b couples the 



gas dynamics of the explosion with the particle acceleration. 
The present form of the solution assumes spherical symmetry. 
In this approximation it is possible to predict the evolution of 
gas density, pressure, and mass velocity, as well as the radii of 
the forward shock and the contact discontinuity, together with 
the energy spectrum and the spatial distribution of CR nuclei 
and electrons, including the properties of their no n-thermal ra- 
diation. A pplied to in dividual young SNRs (see IVolk 1 12004; 



Berezhko. 2005, 2008, for reviews) this theory has successfully 
explained many observed SNR properties. Making use of the 
observed synchrotron emission spectrum from radio to X-ray 
frequencies it permits the derivation of the injection rate of nu- 
clear particles, essentially that of protons, into the acceleration 
process. It has also allowed the determination of the degree of 
magnetic field ampl ification, a process advocated earlier from 
plasma simulations dLucek & Belli bOOOtlBell & Lucekll200ll) . 
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The initia l examples of such a systematic analysi s concerned 



SN 10 06 (IBerezhko et all 120021) and Tycho's SNR dVolk et al 



2002fl It was specifically predicted that magnetic field am- 



plification leads to the concentration of the highest-energy 
electrons and their corresponding synchrotron and IC 7-ray 
emission in a very thin shell ju st behind the forward shock (see 
Fig. 7 of IBerezhko et all 120021) . This concentration, which is a 
result of synchrotron cooling, is observation ally broadened due 



to aprojection effect ( Berezh ko et al. , 2003 1 



As indicated above, the theory has been used in some detail 
to investigate Tycho's SNR (G120. 1+1.4) as the remnant of a 



1 First indications for high magnetic field strengths B = 10 2 — 
10 3 /xG had been obtain ed much earlier for Tycho's SNR by 
iRevnolds & Ellison 1 1 9921 ) from an analysis of the radio spectrum 
alone. 

2 Such filamentary structures have been observed in hard X-rays 
dVink & Lamingll2003l : lLong et alll2003l ; lBamba et alll2003h and are 
nowadays used as a second independent method to infer the mag- 
nitude of the amplified field. As it should be, both methods lead 
to consistent results with i n the errors, whenever both data sets ex- 
ist dVink & Lamingll2003l:lBerezhko & Volkll2004lV61k et al.Ll200l: 
iBalleA 120061 ; IParizot et alll2006h . 
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type la SN in a (roughly) uniform interstellar medium (ISM), in 
order to compare the results with existing data. A stellar ejecta 
mass M e j = 1.4M Q , distance d — 2 .3 kpc, and ISM number 
density N H = 0.5 H-atoms cm" 3 dVolketall 120021 120051) 
were used. For these parameters a total hydrodynamical explo- 
sion energy E sn — 0.27 x 10 51 erg was originally derived to fit 
the observed size R s and expansion speed V s . 

The steep and concave radio synchrotron spectrum shows 
that the forward shock is characteristically modified by 
the pressure of the accelerated nuclear particle compo- 
nent, with the weaker subshock essentially accelerating the 
radio electrons, while higher-ene rgy electrons increasingly 
"see" the overall shock transition ([Reynolds & Ellisonl [1992; 



Berezhko et al., 2002). 



Since the hardening of the particle momentum spectrum be- 
gins at momenta p w m p c, the frequency range, where the cor- 
responding hardening of the spatially integrated synchrotron 
spectrum occurs, brackets the effective mean magnetic field 
strength B,± inside the SNR; here m p denotes the proton mass. 
The magnetic field strength is in addition constrained by the 
requirement that the calculated electron spectrum leads to a 
smooth synchrotron cutoff whose frequency agrees with the X- 
ray observations. From such a comparison of the calculated and 
observed overall synchrotron spectra the field strength resulted 
as £>d ~ 240/iG Assuming the mean magnetic field strength in 
the circumstellar medium of Tycho's SN to be equal to a typi- 
cal ISM- value of 5/iG, and using an overall shock compression 
ratio a < 6 (see below), pure MHD-compression would give 
a much smaller interior field strength < 30^dl A later 
re-analysis, including now independent Chandra observations 
of thin X-ray filam ents as indicators of the outer shock, gave 
B d « 300 ± 60//G dVolk et al.L[2005llah . Overall, together with 
the approximation for the scattering properties of the magnetic 
field fluctuations - summarized in the next section - good con- 
sistency of the predictions of the nonlinear theory with the ex- 
isting observational data was achieved. 

Recently the radius R c of the contact discontinuity (CD), 
separating the shocked ISM and th e ejecta material, has been 
investigated bv lWarren et al. (2005), using Chandra X-ray ob- 
servations. The ratio R c /R s of R c and the radius R s of the 
forward shock (blast wave) is a new physical variable char- 
acterizing Tycho's SNR. It connects the internal dynamics 
with the particle acceleration process. The large mean value 
R c /R s — 0.93 found for this ratio was interpreted as evi- 
dence for efficient CR acceleration which makes the medium 
between those two discontinuities considerably more com- 
pressible. Without accelerated particles, and in a spherically 
symmetric calculation, t he value of the ratio is si gnificantly 
smaller, Rj R a = 0.77 dWang & Chevalieri l200lh . We have 



2007 



2008) that the observed ratio 



demonstrated dVolk et al. 
R c /R s = 0.93 agrees quite well with the results of IVolk et al 
(200l. 



3 Sinc e injection of nu clear particles occurs only at quasi-parallel 
shocks ( I Volk et all 120031) . a pure MHD-compression of a 5fiG up- 
stream field would lead to even considerably lower interior field 
strengths ~ 10/iG in an unmodified shock with a compression ratio 
of 4. This property will be used in subsection 3.1. 



Somewhat later Badenes et al. I d2006l) published a detailed 
comparison between the high-quality X-ray observations from 
XMM-Newton and Chandra of the ejecta emission and cur- 
rent models for Type la explosions. They found that the fun- 
damental properties of the X-ray emission in Tycho are well 
reproduced by a one-dimensional delayed detonation model 
with a kinetic energy E sn = 1.2 x 10 51 erg. Such a value 
is c ompatible with well- known results for deflagration mod- 



els dNomoto et al.L 119841) as well. In any case, this explosion 



energy is rather in the middle of the typical range of type la 
SN explosion energies that vary by a factor of about two (e.g. 



Blinnikovetaill2006t) . In our initial pa- 



iGamezo et all 12004 

per dVolk et all 120071) we had indicated that also such a higher 
kinetic energy of the explosion is consistent with the SNR dy- 
namics. 

Here we shall a priori adopt the value E sn = 1.2 x 10 51 erg 
and perform the calculation of all relevant physical properties 
of Tycho's SNR connected with nonthermal energy production 
and emission in order (i) to find out how well the existing SNR 
data - especially the ratio R c /Rs - are consistent with such an 
explosion energy (ii) to investigate the expected ratio of lep- 
tonic to hadronic 7-ray emission, and (iii) to possibly deter- 
mine the allowed range for the source distance. In section 2 we 
shall outline our approximation for the scattering properties of 
the magnetic field fluctuations, whereas section 3 contains the 
results and the discussion. 



2. Approximation for the scattering properties 
of the system 

Recently iBell d2004l) found that a non-resonant streaming in- 
stability will occur within the precursor of a strong, accel- 
erating shock. The diffusive streaming of accelerated CRs is 
expected to be so strong in this region that a purely grow- 
ing MHD mode appears with a growth rate that is, at least 
at the beginning of the precursor, larger than the growth rate 
of the well-known resonant Alfvenic mode. It is expected that 
due to this non-resonant instability the external magnetic field 
£?ism is amplified within the entire precursor structure, since 
the main effect is produced by the most energetic CRs which 
populate the whole precursor diffusively during their acc elera- 
tion. The Alfven wave excitation in the shock precurso r (Bell 



1978; Blandford & Ostriker, 1978; McKenzie & Volk, 



corresponds to an additional unstable mode that in particular 
leads to a high level of resonantly scattering waves. This mode 
has been argued to dominate the magnetic field energy deep in 
the precursor, after saturation of the non-resonant instability 
(Pell etier et all 120061) . Therefore the overall field amplification 
will be the result of both instabilities operating in the precursor. 

The saturation value B„ r of the non-re sonantly amplifi ed 
magnetic field is given by ( Bell. 12004 Ipellefier et all l2006h 



Bl 



3 V s P c , 

- 1.5 x 10~ 3 P c 



in 2 $ c 



V. 



3000 km/s 



(1) 



where P c is the resulting CR pressure, $ = log(e m ax/TipC 2 ), 
e max is the maximal CR energy, m p is the proton mass, and the 
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value $ = 10 was adopted in the above relation. This value of 
B\ x l (8tt) is by a factor « 3.3 smaller than the value 



B^/(8vr) « 5 x 10" 3 P C . 



(2) 



empirically determined for several SNRs, where the upstream 
(amplified) magnetic field strength Bo in the shock precursor, 
defined as Bq — B&j a, has been derived both from the spa- 
tially integrated radio and X-ray synchrotron spectrum, as well 
as from the thickness of the X-ra y filaments which determ ine 
the downstream field strength B d (Berezhko & Volk, 2006). 

The subsequent growth of the resonant Alfven instability 
towards the precursor takes place in a medium where the effec- 



tive Alfven velocity is given by c\ 



mass density. Following Pell etier et al.1 ( 200 
3.2.2.) one obtains at the subshock 

1/2 

7 



B 



/Airp and p is the 
' (see their section 



(3) 



3000 km/s 

This expression shows that the existing mechanisms are able to 
amplify the magnetic field within the precursor to a strength, 
which considerably exceeds the required value (2) ( Berezhkol 
2008|). In fact, it would be more appropriate to use in Eq.(3) in- 
stead of Bo the value B\ = (a/a s )Bo, which in our model 
is the field strength just ahead of the subshock, where a s 
is the subshock compression ratio. This reduces the numeri- 
cal factor 7 on the r.h.s. of Eq.(3) to about 4 for an object 
like Tycho's SNR, which still substantially exceeds the re- 
quired value. However, this contribution B 2 cs /(8tt) of the reso- 
nant instability to the magnetic field energy density disregards 
any nonlinear wave dissipation. Indeed, within the framework 
of weak plasma turbulence theory, the resonant instability is 
growing freely up to the subshock, where the growth rate goes 

to zero abruptly. 

However, the MHD simulations by iBell (2004), 
cently confirmed with higher numerical resolution 



re- 

Zirak ashvili et al. (2008), show strong dissipation of the ex- 
cited non-resonant instability into the thermal plasma by in- 
duced gas compressions which dissipate in shocks. At the same 
time the field strength is expected to reach a value B nr ^> 
BiSM) discussed ab ove, at which the instability saturates (for a 
different view, see lNiemiec & Pohll 120071) . It is true that these 
simulations do not directly describe the situation in the nonuni- 
form shock precursor with its mixture of non-resonant and res- 
onant fluctuations, where closer to the subshock more and more 
resonant particles appear. Nevertheless the shock dissipation of 
transverse magnetic turbulence into the thermal gas is likely to 
be a generic feature for the mixture of non-resonant and reso- 
nant instability modes, and should in particular also effectively 
dissipate the developing resonant waves. 

This theoretical expectation of strong nonlinear wave dis- 
sipation is consistent with Eq.(2) which shows that the unin- 
hibited growth of the resonant instability leads to a substantial 
overestimate for the magnetic wave energy density at the sub- 
shock. It is therefore a reasonable approximation that at some 
level of magnetic field amplification beyond non-resonant sat- 
uration the nonlinear dissipation by shock formation balances 
the linear growth of the resonant instability, the remaining en- 
ergy going into thermal energy of the precursor plasma. 



The spectrum of CRs produced by a strong, modified shock 
is very hard so that the CRs with the highest energies make 
the largest contribution to the overall CR energy density. These 
most energetic CRs produce a field amplification on their own 
spatial scale which is the precursor size. Therefore CRs with 
lower energies already "see" the amplified field Bo as a mean 
field. Beyond that the strongly excited resonant waves lead 
to strong particle scattering that approaches the Bohm limit, 
where the scattering mean free path locally equals the gyro ra- 
dius in the amplified field B. 

Since the process of magnetic field amplification is not in- 
cluded in our kinetic theory we postulate the existence of a 
far upstream amplified field Bq 3> £>ism that is determined 
from observations of the downstream magnetic field strength 
£?d, and a simplified connection 



B = (p/p d ) B d 



(4) 



between the local field strength B and the local density p. 
Together with the above equation for Bo, this defines the am- 
plified magnetic field everywhere in the shock transition and is 
used in our m odels for the par ticle acceleration and gamma-ray 
emission (e.g. IBerezhkol 120081) . 

A value Bo of the order of 10 £?ism or even l arger, as de 



rived i n the pres e nt pap e r (and for other sources: IVolk et al 



2005; IBerezhkol 120051 120081) . appears quite plausible for 



young objects, in which longer-term dissipation mechanisms 
(Ptu skin & Zirakashvilii 120031 120051) have not yet had suffi- 
cient time to operate. 

The problem is ultimately the scattering (or better, the scat- 
tering strength) of the highest-energy particles at the instanta- 
neous cutoff. This is an intrinsically time-dependent process 
that is not yet resolved i n particle acceleration theor y even for 
unmodified shocks (e.g. Lagage & Cesarskv, 1983), let alone 
for nonlinearly modified shocks jAmato & Blasill2006b . 

This interpretatio n is different from that of 



Ptus kin & Zirakashvilii (120031 120051) who were not yet 



aware o f the non-resonant in stability. It is also differe n t from 
that of lAmato & Blasil d2006l) and Iviadimirov et al.1 d2006l) 
who do not include the non-resonant instability and its field 
amplification effect in their arguments, and rather concentrate 
on the nonline ar effects of the reson ant Alfven wave instability 
in the spirit of Bell & Lucekl (1200 lb . It is finally also different 
from the recent calculation of th e acce leration properties of 
shocks by IZirakashvili & Ptuskin (2008) who base their work 
entirely on the non-resonant instability and its amplification 
properties. 

Our approximation for young objects consists then first 
of all in assuming that the magnetic field strength is ampli- 
fied everywhere around the shock in the sense that all accel- 
erating particles experience the amplified field in the region 
upstream of the subshock. The actual strength of Bo is phe- 
nomenologically determined by iteratively fitting the theoreti- 
cally calculated electron spectrum to the observed synchrotron 
spec trum, and/or by t he analysis of observed X-ray filaments 
(e.g. lVolk et al J, 120051) . The effective downstream field strength 
B d = aB is the result of the MHD-compression of this largely 
perpendicular upstream field. On account of the large ampli- 
tude ( SB 2 ) ~ B 2 of the resonantly scattering waves, the 
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diffusion coefficient is taken as the Bohm limit in the local am- 
plified field B = (p/p,i)Bd, Thi s is clearly a lower lim it to 
the scattering mean free path (e.g. lMarcowifh et"aL , 2006), es- 
pecially for the highest-energy particles whose maximum en- 
ergy is almost certainly overestimated in this way. However the 
Boh m limit seems to be roughly consistent with the observa- 
tions dParizot et al.Ll2006h . 

Secondly, we assume the dissipation rate of magnetic field 
energy into the thermal gas to be equal to the linear growth 
rate caVP c of the resonant Alfven ic wave field, with ca = 
B/yf$mp (see also iBerezhkol 120081) . This does not directly 
include the dissipative limitation of the growth of the non- 
resonantly unstable modes demonstrated by the MHD simu- 
lations and to that extent still underestimates the overall gas 
heating. The contribution of the non-resonant modes to the gas 
heating is only taken into account indirectly through the en- 
hancement of the Alfven velocity in the dissipation term. 

For Tycho's SNR, as for all other cases studied, this gas 
heating lowers the overall compression ratio to values of the 
order of 5 to 7 which seem very plausible numbers in the light 
of the properties of the synchrotron emission . 

It is clear from these arguments that the complexity of the 
magnetic field, even in the asymptotic case of strong shocks, 
can at present only be partially resolved from plasma theory 
alone. The observed synchrotron emission properties of SNRs 
however determine our key parameter Bq which we use to give 
an approximate description as outlined above. This does not yet 
address the questions of the time dependence of this quantity 
(due to the higher shock velocity in the past) or of the late- 
time escape of high-energy particles. In this paper we shall 
assume that Bq is a constant and fixed by the measurements 
at the present time. We note however that even a simplified 
picture with time-independent upstrea m field Bp contains such 
an important element as esc aping CRs (|Berezhko & Krymsky , 
1988tlBerezhko et al.L[l996T) which are the highest-energy CRs 
whose outward diffusion is faster than the SNR shock expan- 
sion. In reality one has to expect this loss process to be even 
much faster due to two factors. First of all, a decreasing B (t) 
due to the shock deceleration makes high-energy CR diffusion 
progressively faster. Secondly, these escaping CR particles that 
exist far upstream from the SN shock produce less efficiently 
magnetic field perturbations and therefore their diffusion co- 
efficient becomes progressively closer to the average diffusion 
coefficient in the ISM. 

3. Results and Discussion 

The dist ance determinations for Tycho's SNR have varied in 
the past. iDe Vaucouleursl (119851) gave a m ost likely range as 



d = 3.2 ± 0.3 kpc. 



Albinson et all dl986b found 2.2+J'g kpc, 
alJ dl995h argued that d could be as large 



whereas ISchwarz et ; 
as 4.5 kpc. Using H Q measurements and assu ming a strong 
gas shock with a compression ratio of a — 4 Kirshner et al. 
dl987l) derived distances of 2.0-2.8 kpc, and lSmithet al]dl99ll) 
extended this to 1.5-3.1 kpc. Very recently iRuiz-Lapuente I 



(2004) estimated d = 2.83 ± 0.79 kpc by comparing Tycho's 
SNR to SN 1006, and d = 2.85 ± 0.4 kpc by using the peak 
luminosity-decline correlation for type la Supernova explo- 
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Fig. 1. (a) Shock (contact discontinuity) radius R s (R c ) and 
shock (contact discontinuity) speed V s (V c ) in units of Rq = 
(4kpc/d) pc and V = 10 4 (4kpc/<i) km/s, (b) total shock (a) 
and subshock (<r s ) compression ratios for Tycho's SNR as func- 
tions of time, calculated for four different distances d. The dot- 
ted vertical line marks the current epoch, t = 435 yrs. The 
observed mean size (filled circle) and speed (fille d square) of 
the sh ock, as determined by radio measurements dTan & Gull . 
1985), are shown as well. 



sions. With the exception of the upper limit of ISchwarz et al 
d 19951) . all these estimates lie below about 3.6 kpc. In such 
a situation we explore here the range of the distances d = 
3.1 — 4.5 kpc. As the lowest value d — 3.1 kpc is taken be- 
cause distances d < 3.1 kpc are inconsistent with existing TeV 
7-ray measurements, as shown below. 

For any given value d we find the density A^h of the ambient 
ISM from a fit to the observed SNR size and expansion speed. 
Due to the relatively small uncertainties of the measurements 
of angular radius R s / d and expansion speed V s /d, there exists 
an almost unique value of A^h for every given distance d, given 
the explosion energy E sn . 

Fig[T] shows the calculations of shock-related hydrody- 
namic quantities for four different source distances in the range 
3.1 < d < 4.5 kpc, together with the azimuthall y averaged ex - 
perimental data from radio observations dTan & Gulll Il985l) . 
The corresponding ambient gas densities, which result from a 
fit to the observed angular SNR size and expansion speed, are 
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log v, 



Fig. 2. Synchrotron emission flux as a function of frequency 
for the same four cases as in FigQ] at the current time. The thin 
dot-dot-dashed line represents the flux expected in the test par- 
ticle limit for an interi or magnetic field s trength = 20 /iG. 
The observed X-ray dAllen et all [1 999) and radio emission 
( Reynolds & Ellison, 1992) are shown as well. 



listed in Table 1 together with the values of other relevant pa- 
rameters at the current epoch. 

The system is just in the transition from free expansion into 
the Sedov phase. Therefore, although inevitably rough, an ana- 
lytical explanation of the resulting relation between Ah and d 
is the following: since in the Sedov phase R s oc (E 8n /Nu) 1 ^ 5 , 
the density Ah oc E sn /d 5 strongly decreases with increasing 
distance d. 



As in our previous studies dVolk et all 1200 2. 2005, 2007) 



we determine the values of the proton injection rate r), electron- 
to-proton ratio K cp and interior magnetic field value from 
a fit to the observed synchrotron spectrum. The softness of the 
observed low-energy radio spectrum - relative to a test particle 
spectrum - requires a proton injection rate r\ — 3 x 10~ 4 in all 
cases for different distances. This implies a significant nonlin- 
ear modification of the shock at the current age of t = 435 yrs 
(see Fig[Tp). 

To fit the measured synchrotron flux a large magnetic field 
Bd ~ 400/iG is required (see Table 1). Only such a high mag- 
netic field strength can provide sufficiently large synchrotron 
losses to fit the observed X-ray flux. As Figf2]shows, one then 
can obtain a fit to the synchrotron radio and X-ray spectra 



which is as good as in our previous studies dVolk et all 120021 
2005), which assumed E sn = 0.27 x 10 51 erg and in lVolk et al 



(2002J) also a slight overestimate of the observed X-ray emis- 



We now compare the theoretical ratio R c /R s of the contact 
discontinuity radius (CD) and the blast wave radius with th e 
observational estimate R c /R s = 0.93 by I Warren et alJ ( 2005b . 
In order to do this, one has to first make the two results com- 



parable. It was argued in IVolk et alJ (2007, 2008) that if one 
starts from a spherically symmetric calculation of the CD ra- 
dius as we do, one has to take into account that the actual CD is 
subject to the Rayleigh-Taylor instability, and thus a correction 
is needed to compare such a 1-D calculation with observations 
of the CD. Here we also correct our calculated CD radius by a 




100 



1000 



t. yr 



Fig. 3. The ratio R c /R s of the contact discontinuity radius to 
the forward shock radius as a function of time, calculated for 
four different source dis tances, together wi th the observation- 
ally determined value of IWarren et alJ 12005). 



factor 1.05. It is clear that the increased compressibility of the 
shocked material can only be due to accelerated nuclear parti- 
cles, because the pressure of the accelerated electrons is always 
negligibly small. 

The comparison of the corrected values R c /R s with the 
experimentally estimated value R c /R s — 0.93 (in Figf3] 
we present that va lue with 2% uncertainties, according to 
Warren et al] d2005l) ) shows quite good agreement in all the 
cases considered (see Fig|3). According to the theoretical cal- 
culation the ratio R c /R s is almost constant in the free ex- 
pansion phase of the SNR evolution, whereas in the adiabatic 
phase the contact discontinuity quickly drops behind the for- 
ward shock and the ratio starts to decrease. According to FigQ] 
and Fig|3] Tycho's SNR is in transit to the adiabatic phase. 
In the case without CR acceleration our corrected value is 
R c /R s ~ 0.89. We note that this value is larger (and the ef- 
fect of CRs corres pondingly smaller) than that assumed by 
Warren et al. d2005l) . because these authors referred to SNR 
parameters acco rding to which Tycho's S NR is closer to the 
adiabatic phase ( Wang & Chevalierll200ll) . 

In order to find a constraint on the distance d and the ISM 
density Ah from 7-ray observations, we compare in Fig|4^ the 
resulting 7-ray spectral energy distribution with the HEGRA 
and Whipple upper limits at TeV energies. The 7-ray fluxes, 
calculated in spherical symmetry, have all been reduced by a 
factor / rc = 0.2, in order to account for the selective injection 
of protons which o perates only over ss 20% of the SNR surface 
dVolketallEool . ] t can be seen that all distances d < 3.3 kpc 
are inconsistent with the HEGRA data^. We also note that at all 
distances considered the expected 7-ray flux is dominated by 
7r°-decay (hadronic) 7-rays (see Table 1 and Figj4j5). However, 
the contribution of these hadronic 7-rays progressively de- 
creases with the increase of the source distance: the hadronic 



7-ray 



flux 



x 



E c N\i/d 2 is proportional to the ISM num- 



4 Given the considerably higher sensitivity of HEGRA, the 
SHALON results contradict the H-CT result. Therefore we shall not 
discuss the claimed detection by SHALON. 
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Table 1. Models Parameters (Tycho SNR) 



d, kpc 


Nh, cm" 3 


(7 


B d , fiG 


B' d> »G 


K cp , 10~ 4 


F*/Fj c 


3.1 


0.5 


5.25 


420 


360 


4.2 


40 


3.3 


0.4 


5.15 


412 


350 


5.1 


27 


4.0 


0.2 


4.85 


388 


320 


9.5 


7.7 


4.5 


0.1 


4.6 


368 


309 


15. 


2.4 




Fig. 4. Spectral energy distributions of the 7-ray emission 
from Tycho's SNR (total spectrum, which is the sum of all 
components, (a), and separately ir° -decay, inverse Compton 
(IC) and non-thermal Bremstrahlung (NB) components, (b)), 
as functions of 7-ray energy e 7 , calculated for different 
source distances, toget her with the upper lim its measured 
by t he HEGRA (H-CT; " 



Aharonian et all 1200 ll) and Whipple 



(W; iBuckley et all 1 19981) Cherenkov te lescopes, the HEGRA 
AIROBICC (HA: IPrahl & ProschL 1 1 997l) upper limit, and data 



presented by SHAWN (S; Sin itsvna et aUl2007l) . The thin dot 



dot-dashed line in Fig|4^ represents the IC spectrum that is ex- 
pected to be produced in the test particle case. 



ber density TYh and to the total CR energy content E c (e.g. 
Berezhko & Volkl 1 1997b . In the Sedov phase E c is an almost 
distance-independent fraction of E sn , and Nh oc d~ 5 . This 
gives then oc d~ 7 , since E sn is fixed. Such a dependence is 
roughly consistent with the results of our calculations given in 
FigH 

Fig H]} shows that for d > 4 kpc the hadronic 7-ray flux 
approaches the inverse Compton (IC) flux. The latter is almost 
insensitive to the source distance d. Therefore its value e 7 F 7 w 



10 -2 eV/(cnrs) at 10 8 < e 7 < 10 lci eV represents the lowest 
possible limit. The contribution of non-thermal bremsstrahlung 
(NB) emission is negligibly small. 



3. 1. Simple IC emission model 

The closeness of the predicted hadronic and IC 7-ray fluxes 
might at first sight also be seen to justify the simple consid- 
eration that a situation is possible or even likely in which no 
nuclear particles are accelerated at all, but only electrons. The 
observed upper limit on the 7-ray flux might then be explained 
by the IC 7-ray emission in the cosmic microwave background 
(CMB) radiation field of the same accelerated electrons which 
also produce the observed synchrotron emission. 

A simple estimate of the expected IC TeV-emission can 
be obtained from the relation e~/F* = vS^m-piy/ {B\/%it) 
which is valid for the gamma ray energy e 7 = 1.5 x 
10 10 (/i^/eV)(10/iG/Bd) eV assuming that IC scattering takes 
place in the Thompson regime. Here w p h sw 0.3 eV/cm 3 is the 
CMB energy density. Since according to Figf2]we have vS v « 
100 eV/(cm 2 s) for hu = 1 keV (y « 2 x 10 17 Hz), the expected 
IC 7-ray energy flux at e 7 = 15(10/iG/i3d) TeV produced by 
the same electrons is e 7 F 7 c w 10(10/iG/i?d) 2 eV/(cm 2 s). 
These relations show that the IC 7-ray flux exceeds the 
HEGRA upper limit unless the interior magnetic field is con- 
siderably smaller than 10 pG. 

In order to illustrate the energy spectrum expected in the 
case of inefficient proton acceleration we present in Fig.4a a 
IC 7-ray spectrum, calculated at £?d = 20 pG an with elec- 
tron energy distribution function f c (e) = Ae~ 2 exp(— e/e ma x), 
which corresponds to the test particle solution. The parameters 
A and e max = 13 TeV in this expression were determined by 
fitting the observed synchrotron flux vS v at v = 1 GHz and 
i/ = 3x 10 18 Hz (see Fig|2]and|4). 

As expected from the results of the nonlinear theory, this 
simple test particle theory which considers an inefficient pro- 
duction of nuclear CRs is quite inconsistent with the existing 
upper limit for the 7-ray production. Since the synchrotron 
emission is expected to come from those regions of the shock 
surface where the field is strongly amplified, the shock is a 
quasi-parallel shock. Therefore the ISM magnetic field which 
should have a strength of 5pG or less, should be MHD- 
compressed by a factor that is significantly less than the shock 
compression ratio. As a result an interior field strength of 
about 10/iG is adequate for this type la SN. In FigH] a field 
strength £>d = 20/iG was chosen. It still overpredicts the 7-ray 
flux by an order of magnitude. Only for £>d ~ 70/iG the pre- 
dicted IC 7-ray flux would be as low as the observed HEGRA 
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upper limrE This result shows that the assumption of electron 
acceleration alone in SNRs has no physical basis and should be 
abandoned. 



3.2. Interior magnetic field strength from X-ray 
filaments 

We have also determined the interior magnetic field value 
from the observed linear thickness L of X-ray filaments (see 
Volk et al. , 2005llai for details). These values are given in 



Table 1 as B' d for different distances. Since the size L in- 
creases proportional to d , the m agnetic field strength B' d cx 
L~ 2 / 3 (Be rezhko & VolkL 120041) decreases with increasing d. 



Therefore there is a growing discrepancy in the ratios between 
Bd and B' d : at d = 4.5 kpc B' d w 310/xG which is already con- 
siderably smaller than B^ w 370//G B' d /Ba is already consid- 
erably smaller than at d = 3.1 kpc. Therefore we might be able 
to constrain the source distance also from above, d < 4 kpc. 



4. Conclusions 



The last result led I Volk et al.1 (120071) to the conclusion that one 
could constrain the source distance also from above, d < 4 kpc. 
However, this is not a strong argument, as the detailed data pre- 
sented in Table 1 show. Therefore, regarding the expected 7-ray 
flux, we can only say the following: the increase in the kinetic 
energy of explosion to E sn = 1.2 x 10 51 erg requires a rather 
high source distance of d > 3.3 kpc to explain the HEGRA 
non-detection. On the other hand, the rather low distance esti- 
mates from independent measurements summarized in the be- 
ginning of this section make it for this increased value of E sn 
even more likely that the actual 7-ray flux from Tycho is "only 
slightly" below the HEGRA upper limit. The strong magnetic 
field amplification implies that the 7-ray flux is hadronically 
dominated, at least above photon energies of 1 GeV. 

From the agreement between our theoretical solutions with 
the measurements of the discontinuity radii we conclude that 
there is new evidence for effective acceleration of nuclear par- 
ticles in Tycho's SNR. A future 7-ray detection would there- 
fore make the case for nuclear particle acceleration in SNRs 
incontrovertible. 



Acknowledgements. HJV is indebted to Wolfgang Hillebrandt for dis- 
cussions about SN explosion energies and to Roger Chevalier for a 
discussion on the modification of the contact discontinuity by the 
Rayleigh-Taylor instability. EGB and LTK acknowledge the partial 
support by the Presidium of RAS (program No. 1 6) and by the Russian 
Foundation for Basic Research (grant 07-02-00221) and the hospital- 
ity of the Max-Planck-Institut fur Kernphysik, where part of this work 
was carried out. 



5 We remark that a similar result is found for SN1006, the other pro- 
totypical type la SN observed, where for Bd ~ the analogously 
calcu lated IC flux equals the H.E.S.S. upper limit dAharonian et all 
l2005t) for the 7-ray flux. Any smaller field Bd < 20^iG would over- 
predict the 7-ray flux. 
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